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ABSTRACT 


Fiber optic communications are the future of U.S. Navy shipboard 
communications. They present tremendous bandwidths with no susecptability to electro- 
magnetic interference (EMI) and outstanding signal-to-noise ratios. Current technology 
uses wavelength division multiplexing (WDM) to allow multiple users on a single fiber 
simultaneously. The optical filters necessary to demcodulate the WDM signal can be 
expensive. A less costly alternative could be achieved if Hadamard-Walsh Code Shift 
Keyed (CSK) encoded signals were used. An optical receiver was proposed, using the 
Mach-Zender coupler, capable of demodulating a Hadamard-Walsh CSK encoded signal. 
Building on previous work studying the characteristics of the Mach-Zender coupler, a 
relationship between the probablility of bit error and signal-to-noise ratio (SNR) was 
developed for a single user and a DPSK optical signal. This relationship was then used to 
develop an understanding of the bit-error rate to SNR relationship for a multiple-user 
CSK optical signal. Using the theoretical performance as a guide, a MATLAB model was 


then constructed to investigate the sensitivity of the receiver to non-ideal components. 
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I. INTRODUCTION 


The majority of data transported through modern U.S. Naval platforms is carried 
on copper wire or coaxial cable. As the amount of data being transferred dramatically in- 
creases, the need for a more efficient methods of transport arises. Given the tremendous 
increase in bandwidth of fiber optics, its use to transport data in naval platforms is be- 
coming standard. Not only does fiber optic communication offer a tremendous band- 
width, it offers many other desirable characteristics that include resistance to electro- 
magnetic interference (EMI) and an excellent signal-to-noise ratio that make its use on 
naval platforms extremely attractive. 

To design the information flow in naval vessels of the future, it is important to 
study methods of optical data transmission and explore ways to improve upon and exploit 
the tremendous bandwidth available. In this thesis an alternative to Wave Division Multi- 
plexing (WDM) will be developed and analyzed. This alternative relies upon the use of a 
specialized coupler, the Mach-Zender coupler, to demodulate the multiplexed optical sig- 
nal into individual signal components. 

The Mach-Zender coupler is composed of two 2x2 optical couplers, a delay line 
and a phase shift element. The signal enters the first coupler input port (port 1) as shown 
in Figure 1.1 and is equally divided between the output ports A and B. Signal a, is then 
delayed one bit period (T) and becomes signal a;. Both signals then enter the second 
coupler are combined and then equally divided to the output ports A and B. The signal 
b; is then phase shifted — 7/2 radians giving the outputs of the Mach-Zender coupler of 
ea and b,. 

Utilizing the outputs of the Mach-Zender coupler, a receiver can be constructed to 
demodulate a single optical Differential Phase Shift Keyed (DPSK) signal as shown in 
Figure 1.2. The outputs of the Mach-Zender coupler are followed by a photodetector 


(PD) to convert the optical signal to an electrical current by an amount proportional to the 


Sa 





a, 


Figure 1.1 The Mach-Zender Coupler 


responsivity (9) of the photodiode. This current signal is further amplified by a transim- 
pedence amplifier of gain Zg and input to an integrator over one period (7). The energy 
of the two channels 1s then compared and a symbol (“1” or “0”’) is chosen based on the 


channel with the greatest energy. 
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Figure 1.2 Demodulation of DPSK Signal Using Mach-Zender Coupler 


By combining multiple couplers with varying delays (e.g. T and 27), a Code Shift Keyed 
(CSK) signal employing Hadamard-Walsh coding could be demodulated. This orthogo- 
nally coded signal would allow multiple stations to use a single fiber optic line simulta- 
neously as is currently being done using Wavelength Division Multiplexing (WDM). 
With WDM, however, there is a significant cost in creating very sensitive filters neces- 
sary to demodulate the signal [Ref. 1]. Figure 1.3 shows a system using multiple Mach- 


Zender couplers that 1s capable of decoding four Hadamard-Walsh coded signals. 
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By increasing the number of Mach-Zender coupler stages, this receiver can be con- 
structed to demultiplex increasingly complex Hadamard- Walsh coded signals. 

In Chapter IJ, the theory behind the Mach-Zender coupler will be discussed, and 
from this theory, a mathematical model of the coupler will be developed. This model will 
include both time and frequency domains. Chapter IJ] will use the Mach-Zender mathe- 
matical model of Chapter IJ to develop a single user DPSK encoded optical receiver. 
From this DPSK receiver a probability of bit-error will be determined based on the sig- 
nal-to-noise ratio at the receiver. Chapter IV will then propose a multi-user Code Shift 
Keyed (CSK) optical receiver and develop the bit error rate for the multi-user receiver. 
Chapter V will then compare the theoretical results of the bit error rate versus signal-to- 
noise ratio in the single user receiver with those obtained by a computer model of the re- 
ceiver. Chapter V will continue with the computer model to conduct a sensitivity analysis 
of the single user receiver. Finally, Chapter VI will discuss the conclusion of this thesis 


and present some areas for future research with these optical receivers. 


Il. DEVELOPMENT OF A MATHEMATICAL MODEL FOR THE MACH- 
ZENDER COUPLER 


A. OVERVIEW OF MACH-ZENDER COUPLER 


In the following analysis, a mathematical description of the Mach-Zender coupler, 
as has been previously developed at the Naval Postgraduate School [Ref. 2], will be re- 
viewed. This model will be used to analyze the demodulation of an optical DPSK signal 
using a single Mach-Zender coupler. 

As previously described, the Mach-Zender Coupler consists of two 2x2 optical 
couplers, a delay line of some delay nT where n is an integer and T is the period of the 
signal, and a device to provide a phase shift of — 2/2 radians. Port 1 of the first coupler 
is the input and Port A of the second coupler and the output port of the final phase shift 
element are the outputs. 

The delay of the Mach-Zender coupler can be accomplished by inserting a length 
of fiber equivalent to the desired delay. Heinbaugh [Ref. 2] proposed using a 2x2 coupler 
with a single input and a single output to obtain a 7/2 phase shift. This mathematical 
model will use a -7/2 phase shift. 

The input signal is defined as s(t) in the time domain and S//) 1s its Fourier trans- 
form in the frequency domain. As the review in this chapter will show, the output in the 
time domain is s(t-T)-s(t) for port A of the second coupler and s(t-7)+s(t) for the output 
of the phase shift element. This is as illustrated in Figure 2.1. 


/2((S(t-T)-s(9) 








/2((s(t-T)+s(t) 


Figure 2.1 Signal Throughput of the Mach-Zender Coupler 


B. DEVELOPING A MATHEMATICAL MODEL 


In this section each element of the Mach-Zender coupler will be analyzed mathe- 
matically and an expression relating to that particular element’s effect on s(t) will be de- 


veloped. 
1. Coupler Analysis 


The optical coupler, Figure 2.2, has four ports labeled 1, 2, A and B. Each port 


can be used as an input or an output . 





Figure 2.2 Optical 2x2 Coupler 


Let the frequency domain inputs of the 2x2 coupler of Figure 2.2 be described by 


the following vectors a and b: 


a, b, 
a ~ |b 
ar : b = ; . Cay 
a, b, 
a, D, 


The input vector @ and the output vector 6 are related by the scattering matrix S 
defined as follows [Ref. 3]: 


Ome | | 

«|| saan 
i iii Qeead (22) 

il “Oia 

where 

b = Sa (2.3) 

b, O 1 jf] 4, 

b| 1/0 0 7 1)@ 
Dele 2s lee ex er) 

b, pe 06 | a. 


Equation 2.2 can be simplified if the matrices of Equation 2.4 are partitioned as 


shown below: 


b, Q 20 eal, 
b| ,/0 0 : Ff 1/4, 
ee | a lee ee tee (2.5) 
b, J2 le J : a, 


Aaa |e I- (a, + ja,)/ V2 a 
See Un | Gre ee yD (2.6a) 
b, ee. J\ 4, _ CEE 

He |, pe alas (2.6b) 


Observe that in the ideal 2x2 optical coupler, no part of the input signal is re- 
flected back to the input ports. Specifically, when the only inputs are a, and a,, the out- 
puts are confined to b, and b, as in Equation 2.6b. That is b, and b, consist of only in- 


phase and z/2 phase-shifted components of a, and @,. 


a> b 





Figure 2.3 Single Input 2x2 Coupler 


To apply this result to a practical use of the Mach-Zender coupler, consider Figure 
2.3 in which one input signal is applied to one input port, and the remaining port is termi- 
nated. Equation 2.6b now reduces to become: 


b, a, / 2 a PAD 
H [a A 7 Re fn (2.7) 


We see that the input a, is decreased by a factor of 1/ V2 at ports A and B and the output 
at port B 1s shifted in phase by z/2 radians from the input. 


2. Delay Element Analysis 


As shown in Figure 1.1 the delay element represents a time shift of 7 in the time 
domain which corresponds to multiplication by e’?” in the frequency domain. The input 
and output of the delay element can be modeled in the time domain as in Figure 2.4 with 
the output @ = a,e’°” in the frequency domain. 





Figure 2.4 Time delay element 


3. Phase Shift Element Analysis 


The phase shift of the Mach-Zender Coupler can be accomplished using a 2x2 
coupler with only one input and one output. The result would be a signal shifted by z/2 
radians. For this model, a phase shift of — 2/2 radians was assumed. The phase shift 


can be modeled in the frequency domain as in Figure 2.5. 





Figure 2.5 Phase Shift Element 


C. INPUT SIGNAL ANALYSIS 


Using the preceding development, the function of the Mach-Zender coupler can 
now be completely described. Consider an input signal, s(t), applied to an input port of 
the Mach-Zender coupler. The unused input port (2) is terminated. Let S(f) be the Fourier 
transform of s(t). The signal can be applied to a Mach-Zender coupler as shown in Figure 
2.6. 





Figure 2.6 Input Signal Through Mach-Zender Coupler 


Using Equation 2.7 the output of the first coupler becomes 


b,|_| S(f)/v2 
i : peel ai 


After the delay line the inputs to the second coupler are obtained by applying the 
delay illustrated in Figure 2.4. The inputs to the second coupler are: 


elo 1 gel aa hg aa 
by} | 0 lpssA/v2] | JS)/V2 | 


The outputs of the second coupler, br" and b , are obtained by applying Equation 
2.6b with b! and b, as the inputs. The following result is obtained: 


Mala | Sec amie) | US (fogs a= Sf) 
3 f=. | =| ie es : (2.10) 
by) 27 WN JSf)/v2 | | F/2S(e"" - SP) 
The addition of the —z/2 phase shift gives the final result in the frequency do- 
main as: 


I : } 0 [; (2(S(f)eP*" — Ad al 9 -S(f)) 
bi) [9 ~sLF/AS(NeP"—SPY] 2LSMeP" + SF) 


| (2.11) 


By taking the inverse Fourier transform of Equation 2.11 the time domain result is 


obtained: 
bi) | 1} s@¢@-T)—-s(¢) 
bi(t)| 2\s(t-T)+s | ax 


By multiplying the matrices used to represent the phase shift, delay and two cou- 
plers (Equation 2.13), we can derive a general equation (Equation 2.14(a), 2.14(b)) where 
the output of the Mach-Zender coupler can be readily determined from the inputs @, and 


aap "lel, in ths, ja, 
alo -i2lLy who afvala ia, 7 


ay, 


Phase 2 delay 12 
Shift Coupler Coupler 
By bene DMG ie (2.14) 
Peal iol Se aoa lee 14a 
BCP Ge ty) [ir 


| URGE ree 
ez ]=4| i )+ He ] | (2.14b) 


a (eP™" +1)+4,(e?* —1) 


D. SUMMARY 


In this chapter the optical 2x2 coupler was shown to combine the signals of its in- 
put ports, split the signal, and send the result to the output port such that both signals 
leave the coupler simultaneously. It was determined that if only one input port was used 


then no part of the input signal was reflected back. It was also shown that the outputs 
from the Mach-Zender coupler are (S(f)e??”"7 —S(f))/2 and (S(f)e?™ + S(f))/2 


1] 


in the frequency domain or, equivalently, (s(¢-—7)—-.s(t))/2 and (s(¢—7)+s(t))/2 in 
the time domain. 

In the next chapter a receiver will be proposed using the Mach-Zender coupler 
that is capable of processing a DPSK encoded signal. Using the results of Chapter II, the 
signal will be traced through the receiver and its strength relative to the receiver noise 
will be analyzed. Finally, a probability of bit-error, P}, will be developed based on these 


results. 


il. DPSK DEMODULATION ANALYSIS 


A. INPUT DPSK SIGNAL 


In Chapter I, an optical receiver capable of decoding a DPSK encoded signal was 
proposed (Figure 3.1). Chapter II then developed a mathematical model (Equation 2.12) 
of the Mach-Zender coupler to be used in this optical receiver. With these results, the 
next step is to analyze the input of a DPSK signal and assess the receiver performance 


versus the input signal-to-noise ratio (SNR). 
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Figure 3.1 Optical DPSK Receiver 


A given DPSK signal will be in one of two forms (over a period of 27 where T is 
the period of the input signal) depending on the information bit being transmitted. For a 
bit “1” being sent where 0 <¢ < 27, the signal will be: 


S,(t) = + Ap, (t — T) cos(2af,t + ) ¥ Ap, (t)cos(27f,t +4) . (3.1) 


For a bit “0” being sent where 0 <¢ < 27, the signal will be: 


S)(t) = tAp-(t — T)cos(27f,t + 8,)+ Ap-(t) cos(27f,.t + @) (3.2) 


where fc is a multiple of //T and p,(t—iT) with i an integer is defined as: 


13 


ee oS a 


0 otherwise. 


pp (t-iT) -| (3.3) 


Due to the symmetry of the signal, an analysis of an input | or 0 will produce the 


same results. To simplify the calculations an input of 0 (Equation 3.2) will be considered. 


B. RECEIVER ANALYSIS 


1. Output of Mach-Zender Coupler 


In Chapter II, the output of the Mach-Zender Coupler was determined in Equation 
2.12. Applying Equation 2.12 to the input signal, given the transmission of a 0 (Equation 
3.2), the following output is obtained where T <t < 27: 


a,(t) os 0 
b,(t) PDA cos(27f.t+0,)+ Acos(27f.t+9,) | (3.4) 


Simplifying this expression yields: 


a,(t) ad 0 
b,(t) | [+ Acos(2af,t + @) (3.5) 


which are the outputs of the Mach-Zender coupler given an optical DPSK input signal. 


2. Photodetector Analysis 


The next stage of the receiver is an optical photodetector (PD) which converts the 
optical signal into an electrical current signal. The PD provides an output current equiva- 
lent to the input peak power (P=A2/2) multiplied by the responsivity ({%) plus the shot 
noise, n,(t) [Ref. 4]. 

Taking our input from the Mach-Zender coupler (Equation 3.5) and applying it to 


our PD gives the following signal 1n the receiver over the interval where 7 <1 < 27: 


14 


a,(t) ” 7, {t) 
b,(t)|~ | R(4? /2) +n, (| ©) 


It is important to note at this stage that the shot noise produced by the two PD de- 
vices (7, (t) and 7, (t) ) are independent of one another. 


3. Amplifier Analysis 


The next stage of the receiver consists of a transimpedance amplifier with gain Zo. 
Along with a gain to the input signal, this amplifier will introduce additive white gaussian 


noise (AWGN), 77amp(t), to the signal. In order to simplify the understanding of noise in 
the receiver a summation of total noise will be necessary. 


Both the amplifier noise, 77gmp(t), and the shot noise, 7, (t) and 7, (1), are mod- 


eled as AWGN with zero mean. If the total noise in each individual channel is 77(t) then: 
M(t) = Nhanp (t) +17, (t) (3.7) 
Therefore, the noise content of each channel following the amplifier stage will be 77, (t) 


and 7,(t) with a power spectral density (PSD) of Ng. Again, the noise in each channel 
will be independent of the noise in the other channel. 


Given this simplifying assumption, the output of the amplifier is simply the input 
(Equation 3.6) plus the amplifier noise (7,,,,(¢) ) multiplied by the amplifier gain (29): 


a,(t) Zn, (t) 
5 a . Pe /2+n, wy G8) 


4. Integrator Analysis 


The next step in the analysis is to take the amplifier output (Equation 3.8) 
and run it through the integrator. This will give us: 
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se J1Zon, lat 


b.(t) ee (Ce), 
[CZ RA? /2) + (Zon, (O) lat 
which will simplify to: 
ie 
om [Zon (lat 
* “4 ' G.10) 


i: . | 
ZRA7T/2+ [[Zn, (Oat 
0 


To continue with the analysis of the integrator output the statistical nature of the 
signal at this point will be examined. Signal a,() is defined as follows: 


a,(t) =a, + N,(t) (3.11) 


Fr 
where a, =0 and N,(/) = [2 7, (t)|dt. Since N,(¢) consists of a constant multiplied 
0 


by zero mean AWGN, its mean is also zero. Therefore a, represents the mean and JN, (f) 


represents the noise. 


Similarly if 5,(¢) 1s defined as: 
b(t) = b, + N,(t) (3.12) 
me 1 — 
where 5, = Z,RA?T/2 and N,(t)= |[Z,7,(f]ar, then 5, represents the mean and 
0 


N,(t) represents the noise. 


Next, for a gaussian random variable such as our zero mean AWGN defined by 
Na(t) and Np(t), the variance (a) is defined by: 


= One ZN (3.13) 


Now that the mean and variance of a,(/) and 6,(¢) has been determined, an 


analysis of the decision process to choose the signal sent can be undertaken. 
5. Decision Variable Analysis 


With the mean and variance of the signals to be input to the decision element de- 
termined, an examination as the probability of an error in signal reception can begin. A 
symbol “0” was assumed to be sent to our device initially, therefore an error will have 
occurred if the symbol “1” is chosen as the output. This probability of choosing a “1” 
when a “0” is sent, Pr{1|0}, will be equivalent to the magnitude of the signal in channel 
a being greater than the signal in channel 5 given that a “0” was sent or Pr{a5>b5|0}. 

If a new variable, Y, is defined where: 


Y=b, -a, (3.14) 
then it 1s apparent that: 
Pria, > b,|0} = Pr{Y <0)0}. (3.15) 
From the analysis of the integrator output it was determined that a, and b,are 
both independent gaussian random variables, and therefore Y, as the difference of a, and 
b;, is also a gaussian random variable. Consequently the mean of Y or Y can be ob- 
tained as follows: 


Y¥ =b,-a@,=Z,RA*T/2. (3.16) 


Similarly the variance of Y(a,*) can also be obtained as follows: 


Cp = Gy one — Agi Zaiiqll = 2 Ay Nall. (aly) 


Pr{1|0 or Pr{¥ <0|0} can be defined by the probability density function of Y as 


follows: 
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Pr{¥ <0} = | f,(l0)ay (3.18) 


where the probability density function, /, ( yl0) , for a gaussian random variable 1s defined 
as: 














ail 

0) = | a | 3.19 
ie (y| ) \2x0, ’ 2° ( ) 

Substituting Equation 3.19 into Equation 3.18 gives: 

0 0 as Pe 
I Oo) 
Pry Y < 0\0; = 0)dy = xO eee LY) 3.20 
{ | } J7,01 )dy eee ’ 20° y ( ) 
| | Cm) _ 
If a new variable x is defined such that x = — then dx = — a Using this 
2 Me 
new variable to substitute into Equation 3.20 gives: 
gla aed aaa x? 

Preece 06) exp] — — |(-o, dx) = —== _~=—*||: exp) -— |. 3.20 
rth ae Lod-F|-a4- Fe J om F fe 029 


Equation 3.21 is now in the same format as the QO function, Q(z), where z= 


Sale 


Therefore: 


Pr{Il0} = f 7) = | AT?) 


Z f2N,T 


pa 
of dona e). (3.22) 





C. PROBABILITY OF BIT ERROR 


The probability of bit error, P, , for this receiver can be determined as: 
P, = Pr(1) x Pri{Oll} + Pr(0) x Pr{ 1/0} . (3.23) 


If a “1” or a “O” is equiprobable (i.e., Pr(1) = Pr(0) ) and, if by symmetry Pr{O|1} 
= Pr{1|0} , then substituting Equation 3.22 into Equation 3.23: 


i aa ar 
P, -of te l= ). (3.24) 


Further, if the peak power (P) is defined as P=A2/2 then Equation 3.24 will simplify to: 


iP 
P= ofp | 5 c | (3.25) 


Defining the voltage signal-to-noise ratio (SNR) as: 


ovr =2- at ny - = RP z (3.26) 
one 2N, ——*Y2N, 














then P, becomes: 


P, = O(SNR). (3.27) 


Figure 3.2 shows this dependence, which illustrates the rapid decline in bit error rate as 
the signal-to-noise ratio exceeds 4 dB. 


D. SUMMARY 


In this chapter a mathematical model for a single user DPSK encoded optical re- 


ceiver was developed using the analysis of the Mach-Zender coupler from Chapter II. 


19 


6 4 = 0 2 4 6 8 
SNR(dB) 


Figure 3.2 Pb vs. SNR 


From this model a relationship of bit-error rate versus signal-to-noise ratio was then de- 
veloped. 


In Chapter IV, a mathematical model of a multi-user CSK optical receiver will be 
developed in a very similar manner to the development in this chapter for a DPSK re- 
ceiver. Again as in this chapter, that model will be used to determine the performance 
(1.e., P, vs. SNR) of the CSK multi-user receiver. 
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IV. CSK DEMODULATION ANALYSIS 


A. INTRODUCTION 


The preceding chapter developed a probability of bit error for a DPSK demodu- 
lated signal. A DPSK demodulated signal implies that only one user can send in informa- 
tion in the optical fiber at a time. In order to allow multiple users to utilize the fiber si- 
multaneously, the signal demodulation must be modified. As was introduced in Chapter I, 
a Code Shift Keyed (CSK) signal employing Hadamard-Walsh coding can be used to al- 
low multiple users on a fiber simultaneously. 

The Hadamard-Walsh code is a straightforward method of creating multiple or- 
thogonal signals. The Hadamard-Walsh basic 2x2 matrix 1s: 


| ala 
4.1 
i of 4.0 


Observe that the correlation of the first row vector with the second row vector is zero. 
The matrix of Equation 4.1 can be built upon to create additional orthogonal vectors as 


follows: 


(4.2) 


It is apparent that the individual] elements of the 2x2 Hadamard-Walsh matrix are used to 


multiply the entire 2x2 matrix and Equation 4.2 can be written as: 


(4.3) 


poh poh — poh 
— 
— 
—" 


Pa) 


Again observe that each row vector is orthogonal (i.e., has a zero correlation) to all oth- 
ers. 

The 4x4 matrix of Equation 4.3 consists of four orthogonal vectors and can be ap- 
plied to a system in which two bits or two users’ information must be decided simultane- 
ously. It is this matrix that will be used to design an optical receiver and develop the 
probability of bit error for two users on the same fiber simultaneously. 

The design of the receiver will be the same as that proposed in Chapter | and is 


shown in Figure 4.1. 


B. INPUT CSK HADAMARD-WALSH ENCODED SIGNAL 


A given two-user CSK Hadamard-Walsh encoded signal, S(t), will be comprised 


of the sum of two orthogonal signals, one each from channel A and channel B. Therefore: 
S(t) =S,(t)+S,(). (4.4) 


S(t) will exist over a bit time of 7=47" where 7” is the period of the input signal. Fur- 
ther: 


= 


T 
= 4. 
; (4.5) 


where T was defined in chapter three as the period of a DPSK signal. The possible signals 


from channels A and B where 0 <t < 47’ are: 


S\,(t) = +Ap,(t —3T')cos(27f,t + 0,) + App (t — 2T’)cos(27f,t + 6, ) 
+ Ap, (t —T')cos(27f,t + 0,)+ Ap, (t)cos(27f,t + 8) (4.6) 


S,,(t) = +Ap,(t —3T')cos(2af,t + 0,)- Ap, (t — 2T') cos(2af,t + 8) 
+ Ap,.(t—T’)cos(27f,t + 0,) — Ap; (t)cos(27f,t + )) (4.7) 


Sop (t) = + Ap, (t —3T')cos(2af,.t + 0,) + Ap, (t — 2T') cos(27f,t + 9) 
— Ap,.(t —T')cos(2zf,t + 0,) — Ap; (t)cos(27f,t + G) (4.8) 


2D 


aé l Le 
Jo' 


mM); 


indjno 
q JouuRYyy) 


«ies 


indjno 
Vv jeuuryr) 


JOAIQI9y Jaydnoz JopuszZ-yoeyy adeIS s[dnny |p ansIy 


(2) Aa (2) ny 


add 


Os es 


st 








}SoBIE] 
3so0yo 









LZ = Aeped 








sajdnoz 
Japusz-yorypy 








d 
jouueys 





().,8 (7S 


she 





jSosIe] 
asooyd 





LZ = Aejaq 





Jajdnoz 
Japusz-yoew 





Vv 
jouueyd 


eS Ges 


(3)#S 


(3)"S 


_L = 4eaq 


sajdnoz 
Jopuaz-yoew 





(1) 


J) 


So4(t) = +Ap, (t — 37’) cos(2af,.t +0,) — Ap, (t — 2T’) cos(2af.t + 8) 
— Ap, (t-—T’)cos(27f_.t + 0,) + Ap, (t)cos(27f,.t +, ) (4.9) 


where the S subscript refers to the bit transmitted and the channel it was transmitted from, 
respectively. For example, S),4 (t) refers to the channel A component of S(t) in which a | 
was transmitted for channel A, and Sog (t) refers to the channel B component of S(t) for 
which a zero was transmitted for channel B. The p,.(t—iT") term is defined as it was in 
Chapter III (Equation 3.3). The frequency f, is a multiple of //T 

Again due to the symmetry of the signal an analysis of all possible inputs is not 
necessary, as each will produce the same result. Also, as each signal in channel A 1s or- 
thogonal to all signals in channel B, there will be no interference between channel A and 
channel B signals in S(t). To aid in the development of a probability of bit error, P,, an 
input of a 1 in channel A and a0 in channel B will be considered. This input corresponds 


to Equations 4.6 and 4.7 respectively and the input signal will be: 


S(t)=S14(t)+Soa(t). (4.10) 


C. RECEIVER ANALYSIS 


1. Output of First Stage Mach-Zender Coupler 


In this multistage receiver the first stage is identical to that analyzed in Chapter I 
(Figure 3.1), substituting T° for T. Again as in Chapter IJ the output of the first stage is 
determined by applying Equation 2.12 to the input signal, Equation 4.9. Given a | in 
channel A and a 0 in channel B the following first stage output is obtained where 
ee Aye 


> | - 2Ap,-(t—3T’) + 2Ap,.(t — 2T’)—2Ap,(t-T’) 


5,(t) 7 C08( If .t+0,) +2Ap_(t—3T’)-2Ap,(t-T’) | (4.11) 


A detailed derivation of this output is contained in Appendix A. 
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2. Output of Second Stage Mach-Zender Couplers 


The second-stage couplers are designed with a 27” delay variable and conse- 
quently the analysis of their outputs is different than that of the first stage. However, with 
only a slight modification, a variation of Equation 2.12 can be used to describe their out- 


puts. Specifically for a delay of 27° Equation 2.12 can be rewritten as follows: 


ca ml is = 2") ef | (4.12) 
b(t)| 2| s(t-—2T’)+ s(t) 


Applying Equation 4.12 to the input of the second stage, Equation 4.11, the fol- 
lowing output is obtained where 37’ <t < 47”: 


San) 0 
Si,(4)|__ 1] 1/2(4Acos(2zf,t + 0) (4.13) 
S.,(t)| 2|1/2(-4Acos(2af,t + 0,)) 
Sip (t) 0 
which can be simplified to: 
Soa(t) 0 


Si,(t) |_| Acos(2zf,t +0) 
Sia(t)| |—Acos(2af,t+,)) | 
Si@) 0 


(4.14) 


These equations represents the output at the second stage Mach-Zender couplers given an 
optical multiplexed CSK Hadamard-Walsh encoded signal. 


2. Photodetector Analysis 


The remainder of the receiver for each channel is identical to the DPSK receiver 
analyzed in Chapter I; therefore, the analysis will be similar. In the next stage of the re- 


ceiver the optical photodetector (PD) converts the optical signal into an electrical signal. 


IS 


The PD provides an output current equivalent to the input power (P=A’/2) multiplied by 
the responsivity, KX, plus the shot noise, 1,(t). 
Taking our input from the Mach-Zender couplers (Equation 4.13) and applying it 
to our PD gives the following signals in the receiver over the interval 37’ <1 < 47": 
Soa(t) nN, (t) 
Si(t)| | R(A*/2)+9, 
Son(t)| |R(A?/2) +n, (| 
5” (t) Ns, (t) 


(4.15) 


It is important to again note at this stage that the shot noise produced by the all the 


PD devices are independent of one another. 
3. Amplifier Analysis 


The next stage of the receiver, as was the case in Chapter III, consists of a tran- 
simpedance amplifier with gain Zo. Along with a gain to the input signal, this amplifier 
will introduce AWGN, Tamp(t), to the signal. As was the case in the analysis of the DPSK 
receiver the total noise in each channel will be summed. Since both the amplifier noise, 
Namp(t), and the shot noise, 7),(t), are AWGN with zero mean, the total noise in each indi- 
vidual channel 1s 77(t) and : 


T(t) = Namp (4) +7, (t) (4.16) 


The noises, 1oa(t), M(t), Nop(t), and 1:p(t), all have a power spectral density (PSD) of No 
and the noise in each channel is independent of the noise in the other channels. 

Given this simplifying assumption the output of the amplifier 1s simply the input 
(Equation 4.14) plus the amplifier noise, Namp(t), multiplied by the amplifier gain, Zo and 
can be written as: 
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Syl, LZoNoa(t) 

Sit) || Z) (RA? 24+, (0) 
Sig] | Z,(RA?/2+15,(0) | 
Se) Lom, (t) 


(4.17) 


4. Integrator Analysis 


The next step in the analysis is to take the amplifier output (Equation 4.17) and 
run it through the integrator. This will give us: 


[[Zom04(0)] at 


Soa) ‘(iz (RA? /2+7, ,(t))] at 
Sin) Milec 0 Ns 


2S (4.18) 
Son (t) {iz (RA? /2+7,, (t))] at 
oD) oe 0 0B 
4T' 
[[Zom» (0) at 
bo i 
which will simplify to: 
4T’ 
[2.70 (t)] dt 
Sew) 4T' 
simgy| | 2oRAT 24 Jom, Ole 
es = ae (4.19) 
sree t 
0 O) | ZR T'/2+ [Zenon (Ola 
Sig (t) jy 
4T' 
[[Zoms (t)| dt 
af 


To greatly simplify the analysis at this point, notice the similarity of the integrator 
output defined in Equation 4.19 to that of the DPSK receiver integrator output defined in 
Equation 3.10. Specifically, the information plus noise channels in both cases (S;/)'(¢), 


2. 


sts Mths 


So, (t) and b,(t)) are identical as are the noise-only channels (S;, (t),So, (¢),and a,(t)). 


Therefore, the statistical nature of the CSK receiver integrator output can be summarized 


as: 


Sees on ton (4.20) 
Se Sone (4.21) 
Sip = Spe PNG (4.22) 
Se See (4.23) 


where the mean of each channel can be described as S,, and S,, (which are both equal 
to zero) and S,, =S,, = Z,RAT'/2. Further, the variances of each channel have identi- 


cal characteristics, and as in Chapter III, each channel represents a gaussian random vari- 


able. Therefore the variance in each channel can be defined as: 
Gee 2 NET (4.24) 


Now that the mean and variance for each channel has been determined, an analy- 
sis of the decision process can be undertaken. 


5. Decision Variable Analysis 


With the determination of the mean and variance of the signals to be input to the 
decision element, an examination of the probability of an error in signal reception can be- 
gin. Due to the symmetry of the signal the probability of a bit error in channel A, PpfAj, 
is equal to the probability of a bit error in channel B, P,{B}. Therefore an error in channel 
A will be considered and the result will be the same for channel B. If a symbol 1 was 


sent in channel A, an error will occur if the symbol 0 is chosen. This probability of 
choosing a 0 when a | is sent in channel A, Pr, {ll}, will be equivalent to the magnitude 


of the signal So, being larger than the signal S;4. A new variable is defined as follows: 
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TS eS oe (4.25) 


It is apparent that to make an incorrect decision Y must be less than zero and therefore: 
Pr, {Ol }= Pr, {So4 > S,a|lt= Pr, < Olt}. (4.26) 


In Chapter II a probability of bit-error was developed (Equations 3.16 to 3.22) 
relating one information plus noise channel to one noise-only channel. Similarly, in this 
CSK receiver, channel A is almost identical to the result derived in Chapter I, with the 
only difference being that the period in channel A is now TJ’ vice T. Using the result ob- 
tained in Equation 3.24 and substituting T’ for T gives: 


-giiga ra 
Pr, {Ol }= 5 74 } (4.27) 


Now using the definition of 7, Equation 4.5, and rewriting in terms of T: 


mili. | iy 
Pr, {Ol f= of 5s a } (4.28) 





D. PROBABILITY OF BIT ERROR 


The probability of bit error for this receiver can be determined as: 


P, =P, {A}= P{B}. (4.29) 


Since the probability of an error in channel A equals the probability of an error in channel 
B, P, can be generalized by considering channel A only. In channel A: 


P, {A}= Pr, (1)xPr, {Olt }+ Pr, (0) Pr, {llOf. (4.30) 
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If a 1 or aO in channel A 1s equiprobable (1.e., Pra(])=Pra(0)=//2) and if by symmetry 
Pra{Oll }=Pra{/\0} then by substituting Equation 4.28 into Equation 4.30: 


P =P {At= o 0 a | (4.31) 


If the power (P) is defined as P=A’/2 then Equation 4.31 will simplify to: 








rad 3 = | (4.32) 


and if the signal-to-noise ratio (SNR) is defined as it was in Equation 3.25, then P, be- 


comes: 


(4.33) 


SNR 
no) 


This result implies a 3 dB degradation in performance when compared to the result for 
DPSK, Equation 3.26. Figure 4.2 shows a plot of the theoretical performance of the 
DPSK and the CSK receivers and at bit-error rates below about 3x10 the DPSK receiver 
is degraded by 3 dB. 


E. SUMMARY 


This chapter built upon the development of the single-user DPSK receiver from 
Chapter II and has shown that the performance of a two-user CSK optical receiver is 3 
dB less than that of a single user DPSK optical receiver. The next step in the study in the 
use of the Mach-Zender coupler in an optical receiver will be to create a computer model 
of the simple DPSK receiver and obtain performance data from the model. If the model 
performance data closely resembles the theoretical data derived in Chapter II then vari- 


ous components of the model can be altered to conduct a sensitivity analysis of the optical 
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SNR(dB) 


Figure 4.2 P, vs. SNR (theoretical DPSK and CSK receivers) 


DPSK receiver. Chapter V will discuss the construction of the computer model and show 


the results of the sensitivity analysis. 


oi 





V. SIMULATION AND SENSITIVITY ANALYSIS 
A. SIMULATION 


1. Design Approach 


The third chapter developed the theoretical probability of bit error for a DPSK 
optical signal demodulated using the Mach-Zender coupler. If a simulation of the optical 
and electrical receiver components can be constructed using MATLAB code and the sta- 
tistical outcomes of the simulated receiver at various signal-to-noise ratios accumulated, 
then a comparison of the simulated performance to the theoretical will confirm the accu- 


racy of the model. 


2. Mach-Zender Coupler Simulation Construction 


Previous work at the Naval Postgraduate School [Ref. 2] has extensively studied 
the physical characteristics of the Mach-Zender coupler. Given this understanding of the 
coupler, a modularized model of the coupler was constructed in MATLAB code. Appen- 
dix B contains a copy of the code for each module. Each element of the coupler was mod- 
eled as an individual function with varying operational characteristics. 

The optical 2x2 couplers were designed with a variable input scattering matrix. 
This allows for varying amounts of coupling of the input signal(s) into the output signals. 
In addition, a delay element was constructed capable of creating any desired delay onto an 
input signal. 

The phase shift element was modeled as another 2x2 coupler as was proposed in 
the phase shift element analysis of Chapter I. A drawback of using another 2x2 coupler 
to phase shift the signal is the 3 dB attenuation of the signal by the coupler. In order to 
balance the final output power of both channels (b/ and b,’ of Figure 2.6) an.attenuation 


element was inserted in the non-phase shifted output. As with the other modules the 


amount of attenuation was a controlled variable. 
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3. Electrical Receiver Element Simulation 


The next stage of the proposed DPSK optical signal receiver converts the de- 
modulated optical output of the Mach-Zender coupler into an electrical signal for proc- 
essing to obtain a digital output signal as described in Chapter HI. Each electrical compo- 
nent was again designed as a module with varying inputs. 

The photodetector can apply the responsivity (‘Xt ) and add additive white gaussian 
noise to its output. The transimpendance amplifier uniformly amplifies the signal and also 
adds additive white gaussian noise. The integrator integrates the signal over one period. 
The relative strengths of the signal in channel A and channel B are then compared to 


make a decision of a “1” or a “Q”. 
4. Statistical Analysis 


The MATLAB model of the receiver was repeatedly run at various SNR to obtain 
the probability of bit error for specific SNR. The simulated results and the theoretical re- 
sults are shown in Figure 5.1. The simulated results differed by about 4/10 of a dB from 
the theoretical. This indicates that the MATLAB model was accurate. 

With the accuracy of the model confirmed, the model can be used to vary the op- 
eration of various components and investigate the sensitivity of receiver performance to 
less than ideal components. Specifically, the amount of delay and the phase shift were 


altered to non-ideal values and new receiver probability of bit errors were compiled. 


B. SENSITIVITY ANALYSIS 


1. Delay Error 


In building the Mach-Zender coupler the delay element was implemented by in- 
serting an extra length of fiber equal in length to one bit period. At 100 Mbps Heinbaugh 
[Ref. 2] showed that tolerances of +0.5% were achievable. Inserting a delay error of of 
0.5% into the MATLAB model showed no difference in P, vs. SNR. In order to 
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Pb vs. SNR (ideal receiver) 


<~—simulated 
theoretical—> 


P = variable 


R=1 
T = 1e-06 
No = 1 





SNR(dB) 


Figure 5.1 P, vs. SNR (ideal receiver) 


acquire measurable data delay, errors of 5%, 10% and 20% were introduced into the 
model. The results are shown in Figure 5.2. With 5% delay error the required SNR for a 
specific Pb increased one half of a dB. A 10% delay error resulted in over a 3 dB degra- 
dation in performance. At 20% delay error the receiver remained a random number gen- 
erator (i.e, Pp=0.5) up through a 13 dB SNR. 
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+— 0% delay error 


Oo — 5% delay error 


x — 10% delay error 


* — 20% delay error 





SNR(dGB) 


Figure 5.2 P, vs. SNR (simulated delay error) 


1. Phase Shift Error 


The phase shift error is defined by the ability of the 2x2 coupler to change the 
phase of the input signal while keeping attenuation as close to 3 dB as possible. In the 
simulation, attenuation was maintained at 3 dB but the effects of phase delay errors of 


10° and 30° were investigated. The results are shown in Figure 5.3. 
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+ — O degree phase error 


O — 10 degree phase error 


x — 30 degree phase error 





SNR(dB) 


Figure 5.3 Pp vs. SNR (phase shift error) 


For a 10° phase shift error receiver performance degradation was minimal and almost 
zero. For a 30° error the receiver degradation was about 2 dB. 


C. SUMMARY 


This chapter reviewed the construction of a computer model in MATLAB code to 
simulate the DPSK optical receiver. The performance of the model was compared to that 
predicted by the theoretical results obtained in Chapter [II and the model performance 


closely resembled theoretical performance. Having established the accuracy of the model, 
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the model was then used to determine the sensitivity of the DPSK optical receiver given 
errors in delay and phase shift elements. Chapter VI will summarize the modeled results 
and discuss areas for future reasearch regarding the use of the Mach-Zender coupler in an 


optical receiver. 
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VI. CONCLUSION 


A. DISCUSSION OF RESULTS 


This thesis proposed the use of the Mach-Zender coupler in a single-user optical 
receiver and mathematically illustrated its theoretical implementation. Next, building on 
this implementation, the theory supporting the use of the Mach-Zender coupler to de- 
modulate multiplexed fiber-optic signals was developed and also illustrated mathemati- 
cally. From this analysis, we created a computer simulation of the single-user optical re- 
ceiver. Comparing the performance of the simulation with the expected theoretical per- 
formance established the accuracy of the model. With an accurate computer simulation 
developed, the characteristics of the model were then altered to investigate the sensitivi- 
ties of the optical receiver to non-ideal components. 

Chapter IV showed theoretically how the optical receiver could be expanded to 
demodulate a two-user multiplexed fiber optic signal. A two-user optical receiver experi- 
enced a 3-dB degradation in performance as compared to the single-user optical receiver 
proposed in Chapter III. The receiver of Chapter IV could be extended to increase the 
number of users multiplexed on a single optic-fiber by extending the number of Mach- 
Zender coupler stages with varying delays. 

Chapter V discussed the receiver simulation using a MATLAB model, and its per- 
formance closely followed theory. This indicated that the model was accurate and that 
adjusting parameters of the model to more closely resemble an actual Mach-Zender cou- 
pler should give reasonable results in the sensitivity analysis. 

The sensitivity analysis conducted with the model showed only one-half dB deg- 
radation from a 5% delay error. As was discussed in the Chapter V, a 5% error was well 
above the 0.5% error achieved in the construction of the Mach-Zender coupler by Hein- 
baugh [Ref. 2]. Therefore, the effects of delay error due to construction on receiver per- 
formance should be minimal. 

The sensitivity analysis then varied the phase shift error up to 30° and a degrada- 
tion of less than 2 dB was experienced. A phase shift error of 30° is large; larger phase 
shift errors are unlikely. Therefore, this receiver design does not appear to be sensitive to 
probable errors in phase shift. 

This sensitivity analysis, conducted using the MATLAB model, indicates that this 


is a realizable receiver with a strong potential of simultaneously demodulating multiple 


5h) 


signals in one fiber optic line. If this method proves to be feasible, it could provide a cost- 
effective alternative to demodulating multiplexed fiber-optic signals over the current 


technique using Wavelength Division Multiplexing. 
B. FUTURE RESEARCH 


Future research on the feasibility of using the Mach-Zender coupler should focus 


On: 


1. Using the existing Mach-Zender couplers to build a DPSK optical-signal re- 
ceiver in order to confirm the sensitivity analysis already conducted with the 
MATLAB model. 

2. Constructing a CSK optical-signal receiver to investigate the reality of de- 


modulating information from two or more users simultaneously on one fiber. 
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APPENDIX A: ANALYSIS OF OUTPUT OF TWO MACH-ZENDER 


COUPLERS FOR CSK ENCODED SIGNAL 


1. SIGNAL INPUTS 


The following two signals are input to the receiver of Figure 4.1 as S(t): 


S,,(t) =+Ap, (t-3T’)cos(2mf,t +0,)— Ap, (t — 2T’) cos(27f,t + 45) 
+ Ap, (t —T’)cos(27f,t +0,)— Ap, (t)cos(27f,t + 4, ) (A.1) 


Sop (t) =+Ap,-(t —3T’) cos(27f,t +0,) + Ap, (t — 2T’) cos(27f,t + 85) 
— Ap, (t -T’)cos(27f,t +0,) — Ap, (t)cos(27f,t + 8,) (A.2) 


due to the orthogonality of the signals and superposition each signal can be 


considered separately and the end results summed. 


2. CHANNEL A INPUT ANALYSIS 


First consider the signal for channel A. To aid in analysis, cos(27f,t+0,) will be 


factored out and left out of the equation until the end result is obtained. The output of 
the first coupler with a delay of T’ for the channel A input will be described as $4 (t) 


and S;(t). The superscript refers to the contributing input channel and the subscript 
refers to the output channel. Therefore from T’ <t <4T’ the output is: 

I 3 
S4(t)= = [Sra(t-T’) — Sya(t)] = 


ol+Apr (t- 47") — Ap, (t —3T’) + Ap, (t — 2T’) — Ap, (t-T’) 
—[+Ap,(t-3T’) — Ap, (t-—2T’) + App (t-T’)— App (OI) 
] , , / 
= 5% 2x [App (t- 30) + App (t- 27) — Apr (t- TY] (A.3) 


4] 


Sf) = S(S140-T) + S140) = 
[+ Ap, (¢ —4T') — Ap,.(t —3T') + Ap,.(t —2T’) — Ap,.(t -T’) 
+ [+ Ap, (1-30) = App (t-2") + App (tT) — 4p (ON 
= =x2x[0]=0 (A.4) 


Equation A.4 is zero and can be disregarded. Therefore, applying the result of 
Equation A.3 to the input of the second stage coupler with a delay of 27” gives an output 
from 37' <t <4T" of: 


Siu" = 5 SAU-T) - SA) = 
o[-Apr 0-51") + Appelt AT") ~ App 0-37") 
=|-Ap-@—38 )+ Ap; =27 )— Ap-@ Sol 
= =x2x(0]=0 (A.5) 


S(O => [S4-T) + SA] 
5 (-ADp (t~ 51") + Apr (C~4")- App t= 37") | 
+[-Apy (t-3T") + Ap, (t- 27") - App (t-T') 
= > x 2x[-Ap..(t —3T’)] (A.6) 


Equations A.5 and A.6 represent the output of the second stage coupler for channel A. 
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3. CHANNEL B INPUT ANALYSIS 


Next consider the signal for channel B. To aid in analysis, cos(27f,t+6,) will be 


factored out and left out of the equation until the end result is obtained. Therefore at the 


output of the first coupler with a delay of T’ for the channel B input from T’ <t < 47” is: 


s°(t) = > [Son(t-T’) — Soa(t)] = 


s[+Ap, (t-—4T’) + Ap, (t —3T’) — Ap,-(t — 2T’) — Ap, (t-T’) 
Sap ee) + Ap. —2 )—Aprir— i )— Apia 
= 5x 2x[-Ap,(t —2T’)] (A.7) 


5®(1) = > (Sia(t-T) + Sa(t)] = 


s[+Apy (t —4T") + Ap,.(t-—3T’)— Ap, (t—2T’)- Ap,(t-T’) 
+[+Ap,(t -—3T’) + Ap, (t —2T’) -— Ap, (t-T’)— Ap, (t)]] 
= =x Do aa Ne ete (A.8) 


where the superscript refers to the contributing input channel and the subscript refers to 
the output channel. Equation A.7 has only one term and that term and a 27” delayed 
version of that term will be out of range of 37°<t<4T’ and can be disregarded. 
Equation A.8 can be applied to the input of the second stage coupler with a delay of 2T’ 
gives an output from 37’ <t< 4T’ of: 


Séq'(0) = 5 (Sal-T?)- So(t)] = 


ol+4, (t- ST’) - Ap, (t-3T’) 
—[+App(t- 37’) — App (t-T)I] 
= 5x 2xI-Ap, (t-37") vex) 
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Siz () = ; [Spt-T) + Spi) = 
s[+ Apr ai) Ape 31) 


pave or )— Apt —7 )) 


= =x 2x[0]=0. (A.10) 


Equations A.9 and A.10 represent the output of the second stage coupler for channel B. 


4. OUTPUTS 


The output of the first coupler in equation 4.11 is as follows: 


sie SiOD+S,0 (A.11) 
S| [sé+s20 | 


and consequently the output of the second coupler in equation 4.14 is as follows: 


SON Sane) 
SO) eleare 
| TSG) 
SCM Sec) 


(A.12) 


Substituting Equations A.5, A.6, A.9, A.10 into equation A.12 will give the result listed 
in Equation 4.14. 
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APPENDIX B: MATLAB SIMULATION PROGRAMS 


This appendix contains the MATLAB code for the various modules of the DPSK 
optical receiver model. 


1. OPTICAL DPSK RECEIVER 


% set input paramters 


f =1e6; % frequency 
lie = 1/f; % signal period 
ie = 1: % length of delay 


on) = 11,53), 1]; % scattering matrix 


ah % attenuation (dB) 
a % responsivity 

No =1; % noise PSD 

Zo = 1; % amplifier gain 


PWR . [8000]; % power levels 
t = linspace(0,T); % time intervals for one signal 


% load number of runs 
n=2000000; 


% query users for number of runs 
%n = input(‘enter number of simulations to run>> '); 


for p=1:length(PWR) 

A = sqr(PWRi(p)); 

% initialize data matrices 

xmtbit = zeros(1,n+1); 

rcvbit = zeros(1,n); 

% initial synchronization signal to receiver 
S1=A*cos(2*pi*f.*t); 

phase = 1; 

% create 1st data signal 
[S2,xmtbit(1)] = onechxmt(f, phase); 
S2 = A*S2; 


% update phase setting 
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if xmtbit(1) == 
phase = -1*phase; 
end 


fort=1:n 
% create 1st data signal 
[S3,xmtbit(i+1)] = onechxmt(f,phase); 
Soi eo, 


% update phase 


if xmtbit(i+71) == 
phase = -1*phase; 
end 


% create next input signal 
Sie loo.oe on), 


% input S to Mach-Zender coupler 
[xa,xb] = mzcoup(S, T, Td,SM,L); 


% input MZ output to RCVR 
{(rcvbit(t)] = revr(xa,xb,R,Zo,No,T); 


% progress signal to next period 


Ss) eyrs 
$2 =S3; 
end 


% resize xmbit to remove last bit 
xmtbit = xmtbit (1:n); 


% initialize error count 
errcnt = 0; 


% compute total nr of errors 
fori=1:n 

errcnt = errcnt + xor(xmtbit(i) ,rcvbit(i)); 
end 


Pb(p) = (errcnt/n); 
SNR(p) = 10*log10(R*PWR(p)*(sqrt(T/(2*No)))); 


end 
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% procedure to save data to .mat file 

clear P S; 

load ideal2dat; 

P=(P,Pb] 

S=[S,SNR] 

clear L AR No Zo SMt T Tdin Pb SNR errcnt S1 S2 S3 rcvbit xmtbit; 
clear f PWR phase p xa xb; 

save ideal2dat; 

clear 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


2. ONE CHANNEL TRANSMITTER 


function [S,xmtbit] = onechxmt(f, phase) 

% 1 CHANNEL DPSK TRANSMITTER one bit output 

% [S,XMTBIT] = ONECHXMT(F) 

% Creates a random transmit (XMTBIT) bit of one or zero. 
% Based on transmit bit it returns a signal 

%  (S) one period in length that is inphase 

% if XMTBIT is a zero and 180 degrees out of 

% phase if xmtbit is one. 


% decide which bit transmitter will simulate 
% signal to send 


vais % compute period of signal 
t=linspace(0,T); % generate 100 time intervals 
X= rand - 0.5; 


ix <= 0 
xmtbit=0; 
S=phase*cos(2*pi"f.*t); 
% compute S for XMTBIT =0 
else 
xmtbit=1; 
S=-1*phase*cos(2*pi*f.*t); 
% compute S for XMTBIT = 1 
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end 


% Naval Postgraduate School, ECE 
% Created by: LT Paul W. Wehner 


3. MACH-ZENDER COUPLER 


function [a,b] = mzcoupsm(S,T, Td,SM1,SM2,SM3,L) 
% MACH ZENDER COUPLER with varied scattering matrices output 
% For agiven input signal, S (with a period T, 

% three periods in length), returns two signals, 

% A andB (both one period in length) where 

% Az=S(t-T_actual)+S(t) and B=S(t-T_actual)-S(t). 
% T_ACTUAL indicates the actual length 

%  (i.e.,seconds) of the coupler delay 

% element and SM indicates the 

% scattering matrix of the individual #'rd 

% 2x2 optical couplers. L is the attenuation 

%  inchannel a to equalizer signal strength 

% lost in channel b due to the last phase 

% shift. , 


determine # discrete quantities in 
a signal period (T) of S 
n = length (S)/3; 


oO 
o- o- 


% demux signal S in 1st 2x2 coupler 
[a1,b1] = optcoup2(SM1,S); 


% delay signal in channel a 
[a2,d] = delay(a1,T,Td); 

% adjust length of channel a and 
% Channel b signals 

a3 = a2(d-n+1:d); 

b2 = b1(d-n+1:d); 


% 2nd 2x2 coupler output 
[a4,b3] = optcoup2(SM2,a3,b2); 


% phase shft element output 
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[q,b4] = optcoup2(SM3,b3); 


% attenuate channel a 
[a5] = atten(a4,L); 


% define MZ coupler output 
a = real(a5); 
b = real(b4); 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


4. OPTICAL COUPLER 


function [a,b] = optcoup2(SM,S1,S2) 

% optical 2x2 coupler 

% OPTCOUP2(SM,S1,S2) 

% For agiven input signal, S1, and a possible 
% second signal S2, provides two outputs, 

% AandB. Both outputs are attenuated by a 
% factor of 3dB and multiplied by the scattering 
% matrix, SM. SM must be [2x2]. 


% determime if second signal S2 was input 
if nargin == 
$2=zeros(1,length(S1)); 
elseif length(S1) ~= length(S2) 
disp(‘ERROR - Input arrays (S1 and S2) must’) 
disp(’ be of the same length’) 
a=0; 
b=0; 
return; 
end 


% define input matrix 
= ($1;52]; 


% multiply by scattering matrix 
X = SM*y; 


% attenuate output by 3dB 
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Xe Sanic): 


% define individual outputs 
ee), 
b = x(2,:); 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


5. DELAY ELEMENT 


function [a,d] = delay (a,T, Td) 

% DELAY ELEMENT 

% DELAY(A,T,Td), fora given signal A (A must 
% be three periods in length with a period of 

% 1 seconds) outputs a new signal A delayed Td 
% seconds. Indicates at which element 

% Value, N, in a the delayed vector begins 

% where the vector Is listed right to left 


% Create working copy of a 
b= 4; 


% determine # discrete quantities of one 
% periodofa 
| = length(a)/3; 


% compute delay as percentage of actual delay 
perc_Td = Td/T; 


% equate percent delay to # of discrete 

% quantities of a 

j= round(i* perc_Td); % quantized delay 
K = (3%) - j; % new start of a 


% zero out delayed discrete elements of a 
for | = (3*1):-1:(k+1) 

a(l)=0; 
end 
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i=" 


% overwrite a with delayed version of itself 
for n= k:-1:(2*i-j+1) 

a(n)=b(m); 

m=m-1; 
end 


d=k; 


% Naval Postgraduate School, ECE 
% Created by: LT Paul W. Wehner 


6. ATTENUATOR 


function [a] = atten(a,L) 

% ATTENUATOR 

% ATTEN(A,L) 

% For given input signal, A, performs 
% attenuation on signal equal to 

% dB. 


a =a/(sqrt(10“(L/10))); 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


7. ELECTRICAL CIRCUIT ELEMENTS 


function [rcvbit] = revr(xa,xb,R,Zo,No,T) 

% OPTICAL DPSK SIGNAL RECEIVER OUTPUT 
% [RCVBIT]=CRCVR(XA,XB,R,ZO,NO,T) 
% Forgiven input signals, XA and XB 

% this applies photodetection 

% proportional to the responsivity (R) to 
% convert the optical signal to electrical, 
% amplifies the electrical signal 

% with again of (Zo) and then 

% integrates the signal over a 

% single time period (T) to obtain signal 


S| 


% energy in each branch. A comparison of 
% two branches is then made to 

% determine the transmitted bit. Noise is 
% added in the photodetector and 

% the amplifier with a power spectral 

% density of No/2. 


% inut optical signal to photodetector 
xa1 = pd(xa,R,No, 7); 
xb1 = pd(xb,R,No,T); 


% input signal to transimpedance amplifier 
Xa2 = tranamp(xa1,Zo,No, T); 
xb2 = tranamp(xb1,Zo,No, T); 


“% integrate both channels over T 
xa3 = intg(xa2,T); 
xb3 = intg(xb2, T); 


% compare channels to decide received bit 
if xa3 > xb3 


revbit = 1; 
else 

rcvbit = 0; 
end 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


8. PHOTODETECTOR 


function [a] = pd(x,R,No,T) 

% PHOTODETECTOR 

% PD(x,R,No,T) 

% For given input signal, X, performs 

% photdetection on signal power proportional 
% tothe responsivity, R, in amps/watt. 

% Assume signal is output over a 1 ohm 

% load such that amps correlate directly to 
% volts. Adds photodetection shot noise to 
% Output signal A with PSD(No/2) as defined 
% from user input. 
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% Create shot noise 
s = randn(1,length(x)); 


% adjust shot noise PSD 
amp = sart(No/T); 
S=amp*s; 


% compute signal power 
Px = (x.42)./2; 


% multiply input power by responsivity 
% to get output in amps 
Ax = R * Px; 


% add shot noise to create output signal 
a=Ax+sS; 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


9. TRANSIMPEDANCE AMPLIFIER 


function [a] = tranamp(x,Zo,No, T) 
% TRANSIMPEDANCE AMPLIFIER 


% 
% 
% 
% 
% 
% 
%o 


TRANAMP(x,Zo,No,T) 

For given input signal, X, performs 
uniform amplification on signal 

with gain Zo. Adds noise to signal prior 
to amplification with PSD No/2. 

Assume signal is output over a 1 ohm 
load such that amps correlate directly to 
volts. Adds photodetection shot noise to 
output signal A with PSD(No) as defined 
from user input. 


Create amp noise 


s = randn(1,length(x)); 


% adjust amp noise 
amp = sqrt(No/T); 
S=amp*s; 
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% add amp noise to input signal 
aA=X+S; 


% adjust output by gain, Zo 
a=ZO* a; 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


10. INTEGRATOR 


function [a] = intg(x, T) 

% INTEGRATOR 

% INTG(x,T) 

% Forgiven input signal, X, performs 
%  trapazoidal integration on signal 
% over one time period T. Assumes 
% elements of X are equally spaced 
% fromOto T. 


a = (T/length(x)) *trapz(x); 


% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 


11. Q FUNCTION 


function [Q,P] = qfn(x) 

% This function computes the Gaussian Q-function 

% using the rational approximation 26.2.17 of 

% Abromowitz and Stegun. For sufficiently large 

% arguments, the asymptotic expansion 26.2.12 is used. 
% The rational approximation is accurate to within 

% 7.5e-8. Q is the value of the Q-function (integral of 

% Gaussian pdf) of zero mean and unit variance from the 
% input quantity to infinity, and P = 1 - Qis the 

% cumlative distribution function. The call is 
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% [Q,P] = qfn(x). 


% R.E. Ziemer 
J 4/17/94 


b(1)=0.3 1938 153; 
b(2)=-0.356563782; 
b(3)=1.78 1477937; 
b(4)=-1.821255978; 
b(5)=1.330274429; 
p=0.23 16419; 
y=abs(x); 
T=1/(1+p*y); 
TT=zeros(1,5); 
for r=1:5 
m= 4r: 
end 
Z=exp(-y*2/2)/sqrt(2*p1); 
if y<4 
=7,*b* IT: 
else 
Q1=(Z/y)*(1-1/y42+3/y*4); 
end 
if x>0 
Q=Q1; 
else 
Q=1-Q1; 
end 
P=1-Q; 


ae, 
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